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Introduction. While several studies beginning in the 
Apollo-era and continuing up to the present have tried to detect 
the lunar core, either geochemically, geophysically or geolog- 
ically [e.g. 1, 2, 3, 4], it has so far proven somewhat elusive. 
The unambiguous detection of the lunar core is of prime im- 
portance as it holds the potential of distinguishing between the 
various theories for the formation of the Moon. The theory 
which currently enjoys the greatest success is the giant impact 
model, which has the Moon forming about 4.5 Gyr ago from 
the debris produced when a Mars-sized proto-planet collided 
with the proto-Earth [e.g. 5, 6], Simulations reveal that the 
material front which the Moon is made up contains very lit- 
tle iron and consequently a lunar core, if it exists, should be 
small. 

In a recent investigation of more than 30 years of lunar 
laser ranging (LLR) data, Wiiliams et al. [6] detected a dis- 
placement of the Moon’s pole of rotation, indicating that dis- 
sipation is acting on the rotation. Two effects were proposed, 
one being due to solid-body tides raised by the Earth and the 
Sun and the other as stemming front a fluid core. Williams 
et al. found that in order to account for all dissipation terms 
a model including a liquid core had to be invoked. Using the 
approximate boundary layer theory of [7], maximum radii of 
352 km for a liquid Fe core (p= 7.0 g/cnt 3 ) and 374 km for a 
Fe-FeS eutectic (p= 5.3 g/cnt 3 ) were found, respectively. This 
argument that tides plus a fluid-core/solid-mantle interaction 
satisfactorily explain the lunar rotational dissipation data, is 
the main point put forward by Williams et al. in arguing for a 
present-day molten lunar core. 

It is the purpose of the present paper to infer information 
about the deep interior of the Moon, not least to try to sub- 
stantiate the existence of a liquid core, by inverting several 
geophysical observations obtained through LLR analysis and 
Doppler tracking of the Lunar Prospector spacecraft, specif- 
ically, the second degree tidal Love number, mass, moment 
of inertia and the tidal quality factor. Or paraphrased, the 
approach taken here, and detailed in [8], is to explore, in a 
quantitative manner, through a rigorous inversion, the impli- 
cations of the solution parameters found by Williams et al. 
pertaining to lunar geophysics (notably fe). 

Method of Analysis. The inverse problem dealt with 
here of obtaining information on the lunar density and 5-wave 
velocity profile front the four numbers follows our earlier in- 
vestigations by employing an inverse Monte Carlo sampling 
method (see [8] for details). The basic question that we want 
to answer here, is the following 

• how likely is it. given the observed data, their uncer- 
tainties as well as prior information, that the Moon has 
a molten, partially molten or solid Fe, Fe-S alloy or 
dense silicate core (the latter may contain elevated Fe 
and Ti abundances)? 


Our model of the Moon is divided into 5 spherically symmetric 
shells. Each of these shells is physically described by the fol- 
lowing set of parameters, density, bulk and shear moduli, local 
dissipation and layer thickness. In wanting to answer the above 
question, certain conditions are imposed on some parameters, 
notably the shear modulus in the innermost layer (taken to be 
the core). As the shear modulus is the parameter that phys- 
ically distinguishes between solid and liquid material, it can 
only vary within certain intervals, which in terms of 5-wave 
velocity translates into vs SE [0; 0.5] km/s (liquid core) and 
vs > 3.5 km/s (solid core). In sampling the density profile, the 
assumption is also made that it increases as a function of depth. 
Apart from this, prior information concerning the remaining 
parameters is in the form of broad homogeneous probability 
density distributions (ppd’s), reflecting the fact that little is as- 
sumed known about the interior of the Moon. We use Bayesian 
hypothesis testing in the form of the Bayes factor to evaluate, 
in the light of the data and the prior information chosen, which 
scenario, i.e. whether the core is solid or liquid, is the most 
likely. Bayesian hypothesis testing concerns itself with having 
to distinguish between competing hypotheses given data and 
prior information, such as is the lunar core molten (hypothe- 
sis i) or solid (hypothesis j)I The two hypotheses obviously 
correspond to different areas of the model space, Mi and Mj , 
respectively. If {mi, m 2 , . . . , m*,} is a sample from Mi 
of L( m) containing k points, and {mi, m 2 , . . . , ihf} is a 
sample from Mj of L(m) containing l points (L( m) are 
integrated likelihoods, i.e. L( m) = J M L(m)dm), then the 
Bayes factor can be estimated using the expression 
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Consequently, the interpretation of the Bayes factor is such 
that if Bn, > 1, hypothesis j is relatively more plausible given 
d; if, on the other hand, B^. < 1, then hypothesis i has in- 
creased in relative plausibility. On a more technical note it 
just means running separate algorithms corresponding to the 
different hypotheses of interest and then using all the samples 
from these runs to evaluate the Bayes factor as given in the 
expression above. 

Results. The Bayes factor was found to be Tfe=0.00014 
and as Bjj <1, it signifies that hypothesis i is more plausible 
than hypothesis j or in words, that given data and prior infor- 
mation, a liquid core is more probable than a solid core. Our 
results are shown below in the form of posterior 2D marginals, 
showing cross correlations between several parameters. These 
figures highlight most probable state, size and density of the 
individual layers making up the Moon. From figure 1 it can be 
seen that there exists a high probability for a core with a radius 
of around 350 km and density ~7.2 g/cm 3 . Other possibilites 
are of course also possible, most notably there seems to be a 
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large probability for a core with a density of ~6 g/cm 3 and 
radius around 400 km, although it has to be noted that the 
peak is much narrower and not as high as the other one. From 
figure 2 depicting the correlation between shear wave velocity 
and core radius, a notable spike is apparent at about r=350 km 
and v s very close to zero. Figure 3 examines the correlation 
between the shear wave velocity and density in the individual 
layers and these indicate with a high probability a core with a 
density around 7.2 g/cm 3 and v s close to zero as consistently 
inferred hitherto. 

In conclusion, the most likely outcome of our study, based 
on the data, their uncertainties and prior information, is a cen- 
tral core with a most probable 5-wave velocity close to 0 km/s, 
density of ~7.2 g/cm 3 and radius of about 350 km. This is 
interpreted as implying the presence of a molten or partially 
molten Fe-core, in line with evidence presented earlier using 
LLR regarding the dissipation within the Moon. 


density of ~6 g/cm 3 . 




Figure 1 : 2D posterior probability density showing the corre- 
lation between density and thickness for each individual layer. 
Layer #1 corresponds to the core and layer #5 to the crust. 
Colours indicate probability, with dark red signifying the most 
probable and dark blue the least probable solution. A signifi- 
cant amount of probability centered at a core radius of ~350 
km and density ~7.2 g/cm 3 is apparent. Another possibility, 
although less probable, is also visible at a radius ~400 km and 


Figure 2: 2D posterior probability distribution showing the 
correlation that exists between the parameters vs and thick- 
ness for each individual layer. Layer #1 corresponds to the 
core and layer #5 to the crust. Colour coding as in figure 1. 
A solution at a radius of ~350 km and v s close to 0 km/s is 
indicated. 



Figure 3: 2D posterior probability distribution showing the 
correlation that exists between the parameters v s and density 
for each individual layer. Layer #1 corresponds to the core 
and layer #5 to the crust. Colour coding as in figure 1. A so- 
lution with a large probability at a density of ~7.2 g/cm 3 and 
v s close to 0 km/s is apparent. 


References. [1] C. T. Russell et al., Proc. Lunar Planet. 
Sci. Conf. XII, 831, 1981. [2] K. Righter, Icarus, 158, 1, 2002. 
[3] L. L. Hood, Geophys. Res. Lett., 26, 2327, 1999. [4] J. G. 
Williams et al., J. Geophys. Res., 106, 27933, 2001. [5] A. G. 
W. Cameron, in Origin of the Earth and Moon, Univ. Arizona 
Press, Tucson, 2000. [6] R. M. Canup, Icarus, 168, 433, 2004. 
[7] C. F. Yoder. Icarus, 117, 250, 1995. [8] A. Khan et al., J. 
Geophys. Res., 109, E09007, doi: 10.1029/2004JE002294. 

















